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Modified kagome physics in the natural spin-1/2 kagome lattice systems — kapellasite 

Cu3Zn(OH)6Cl2 and haydeeite Cu3Mg(OH)6Cl2 
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The recently discovered natural minerals Cu3Zn(OH)6Cl2 and Cu3Mg(OH)6Cl2 are spin 1/2 sys- 
tems with an ideal kagome geometry. Based on electronic structure calculations, we develop a 
realistic model which includes couplings across the kagome hexagons beyond the original kagome 
model that are intrinsic in real kagome materials. Exact diagonalization studies for the derived 
model reveal a strong impact of these couplings on the magnetic ground state. Our predictions 
could be compared to and supplied with neutron scattering, thermodynamic and NMR data. 

PACS numbers: 71.20.Ps, 75.30.Et, 91.60.Pn 



For decades, low-dimensional spin 
broad interest due to their intriguing. 



systems attract 
unusual ground 



states (GS) such as helically ordered, spin Peierls, spin- 
hquid or resonating valence-bond GS's [H, [3, SB- These 
unusual GS's are typically driven by competing inter- 
actions or geometric frustration. Two-dimensional (2D) 
quantum spin systems are of particular interest because 
the competition between quantum fluctuations and inter- 
actions seems to be well balanced, and fine tuning of this 
competition may lead to zero-temperature transitions be- 
tween semi-classical and quantum phases There are 
several examples for strongly frustrated 2D quantum spin 
materials, e.g. PbVOa or SrCu2 (603)2 0, which 
can be well described by a frustrated spin-1/2 Heisen- 
berg model. Such 2D quantum magnets are at present 
most suitable objects for the comparison between theory 
and experiment. 

A simple but very challenging realization of a geometri- 
cally frustrated quantum magnet is the spin-1/2 Hcisen- 
berg antiferromagnet (HAFM) on a kagome lattice. The 
kagome HAFM attracts much interest due to its unusual 
classical and quantum GS's and low-temperature ther- 
modynamics, see e.g. Refs. H, H, [13, El 12, 13, and 
also due to potential applications of a possible quan- 
tum spin-liquid state [l4 . ITsj . The recent discovery of 
a natural spin-1/2 kagome compound Cu3Zn(OH)eCl2 
(mineral herbertsmithite [l6|) a. nd a subsequent synthe- 
sis of good-quality samples [l7| have spurred both ex- 
perimental m, [iil and theoretical [nl, [H, [H inves- 
tigations of this frustrated magnetic system. The ex- 
perimental results were quite unexpected: Curie- Weiss 
behavior with a rather large Q, an upturn in magnetic 
susceptibility at 75 K and no spin gap down to 100 mK 
are far from being consistently described by theory. The 
main obstacle for theoretical studies is the structural 
Cu — Zn disorder within this compound (20j , which ham- 
pers the kagome physics, but encourages the search for 
new materials. A very recent discovery of two isostruc- 
tural spin-1/2 kagome systems — the minerals kapel- 



lasite Cu3Zn(OH)6Cl2 ([21[, a metastable polymorph 
of herbertsmithite) and haydeeite Cu3Mg(OH)6Cl2 |22j . 
widens the range of possible investigations. These sys- 
tems are of great potential interest because (i) no cations 
are located between the planes, thus less coupling be- 
tween kagome layers is expected though the interlayer 
distance is reduced by about 1 A (ii) the presence of 
two isostructural compounds should allow a systematic 
study of additional exchange couplings beyond the orig- 
inal kagome model. We have performed a theoretical 
electronic structure study within density functional the- 
ory (DFT) and estimated the exchange parameters of a 
corresponding Heisenberg model. For this spin model we 
have calculated the classical GS and for a finite lattice of 

= 36 sites the quantum spin-1/2 GS. 

The DFT calculations were performed using a full- 
potential nonorthogonal local-orbital scheme (FPLO ver- 
sion 6.00-24) [2^ within the local density approxima- 
tion (LDA). The Perdew and Wang parameterization 
of the exchange-correlation potential was chosen for the 
scalar relativistic calculations [2J|. The default basis set 



was used. The strong on-site correlations of the Cu d- 
electrons were taken into account using the LSDA-|-fJ 
method [2^. Well converged i-meshes of 124 points for 
the conventional cell and 75 points for the supercell in 
the irreducible wedge were used. 

The hexagonal crystal structure of both minerals con- 
sists of layers (Fig. [T]) perpendicular to the c direction. 
These layers arc built by a kagome lattice of corner- 
sharing CUO4 plaquettcs, which are tilted with respect to 
this plane, and ZnOg (kapellasite) or MgOg (haydeeite) 
octahcdra bridging the "ring" of six CUO4 plaquettes. 
The Cu — O — Cu angle between two neighboring plaque- 
ttes is close to 105°, providing considerable ferromagnetic 
(FM) contributions to the exchange due to the vicinity 
to 90°. The kagome layers are separated by CI atoms, 
which are bonded to H atoms that stick out of the layers. 
The experimentally defined H position for haydeeite [2^ 
yields the unusually short O — H distance of 0.78 A, the H 
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FIG. 1: (Color online) The hexagonal crystal structure of 
kapellasite and haydeeite: CUO4 plaquettes form a buckled 
kagome layer bridged by ZnOe/MgOe octahedra. The inter- 
layer space is filled with CI and H (not shown) atoms. 
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0-H distance (Angstroem) 

FIG. 2: (Color online) Total energy per H atom given by LDA 
calculations for kapellasite (Zn) and haydeeite (Mg). The zero 
energy corresponds to the calculated equilibrium distance and 
is marked with a gray line. The experimental value of the O — 
H distance in haydeeite is shown with a dotted line. Inset: 
exchange Ji from the supercell LSDA+17 calculations as a 
function of the O — H distance. 



position in kapellasite has not been reported. To account 
for this structural peculiarity, the H position was relaxed 
with respect to the total energy, we show below that it 
has a dramatic impact on the exchange. Throughout the 
paper, we use the optimized H position [27[ yielding ~1.0 
A for the O— H distance (Fig.H]). 

Our LDA calculations yield a valence band with a total 
width of 6 — 7 eV for both compounds with three bands 
crossing the Fermi level £f according to the three Cu 
atoms per unit cell (Fig. [3]) . The valence bands of hay- 
deeite and kapellasite have two pronounced differences: 

(1) the rather localized d-states of Zn (between —6.5 and 
—4 eV) contribute to the valence band of kapellasite while 
Mg states have negligible contribution for haydeeite and 

(2) the width of the separated band complex at ep is 
slightly different. Nevertheless, the same model can be 
applied for the description of low energy excitations. 

The band structure (Fig. [3]) reveals that the dispersion 
perpendicular to the kagome planes (along F — A) is very 
small, pointing to a pronounced 2D character of the sys- 
tems in accordance with our expectations. The presence 
of states at Fermi level yields a metallic GS, contrary to 
the insulating behavior typical for undoped cuprates [23] • 
This discrepancy originates from the strong on-site cor- 
relations of the Cu 3d electrons, insufficiently described 



FIG. 3: (Color online) TB model (squares) fitted to the band 
structures (solid lines) of kapellasite (a) and haydeeite (b). 



by LDA, and can be accounted for by adding the miss- 
ing Coulomb repulsion in a model Hamiltonian or in the 
LSDA-f- U approximation. Though LDA fails to describe 
the correlations correctly, it is known to provide reliable 
values of transfer integrals [2^ which can be used for a 
model analysis of the magnetic excitations. To define 
the relevant orbitals for the low energy excitations we 
analyzed the density of states (DOS) by applying local 
coordinate systems for all orientations of CUO4 plaque- 
ttes and calculated the local DOS and band weights. The 
analysis revealed that the bands at Fermi level belong to 
local Cu 3da;2_j,2 and O 2p^ orbitals, i.e. the standard 
cuprate scenario with a half-filled antibonding dpa* band 
is realized. Therefore, an effective one-band model, al- 
ready applied for similar materials SEI, is appropriate 
to describe the magnetic excitations in these systems. 

Three dpa* bands per unit cell lead to a 3x3 matrix 
representing the tight-binding (TB) Hamiltonian. The 
number of transfer integrals, included into the TB model 
was picked to get a good fit of the LDA bands which 
could not be considerably improved by inclusion of fur- 
ther parameters. The fits shown in Fig. [3] were achieved 
using ten transfer integrals, though only four of them 
(Fig. [1]) were larger than 10 meV. To check the stability 
of the leading terms, we subsequently decreased the num- 
ber of parameters in our model. Basing on these results, 
we estimate less than 10% uncertainty in our values for 
the leading four transfer integrals depending on the cho- 
sen TB Hamiltonian. Thus, we can restrict ourselves to 
analysis of the leading terms. In order to estimate the 
antiferromagnetic (AF) exchange, the transfer integrals 
were mapped to an extended Hubbard model and subse- 
quently to a Heisenberg model with jf^ = Atf/Ueff (s^ . 

The leading AF exchange in both systems is the 
nearest- neighbor (NN) exchange Jf^, the second largest 



is the exchange along diagonals of a kagome lattice, J^^ 
(see Fig.[T]). The relevance of the latter exchange is rather 
unexpected: while the pure kagome model includes Ji 
only, its modifications usually contain only the second 
neighbor J2 [s^. In our case, we find J^-'^ (and also J^^f) 
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FIG. 4: (Color online) Exchange integrals of kapellasite and 
liaydeeite as a function of Coulomb repulsion U. 
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TABLE L Transfer and exchange integrals of kapellasite and 
haydeeite. All values are given in meV. The values of transfer 
integrals are taken from the TB model. The AF exchange 
is calculated via subsequent mapping the transfer integrals 
to the extended Hubbard and Heisenberg models. The total 
exchange is taken from LSDA+ U total energy calculations of 
supercells. J^'^' is the difference between J and J^^ . 



smaller than 0.5 meV for both systems, thus these terms 
can be neglected in the following discussion. The inter- 
plane coupling is much smaller than 0.1 mcV, showing 
that the systems arc almost perfect 2D magnets. There- 
fore, a Ji-J d model should be appropriate to describe the 
magnetism in good approximation. 

The values of the total exchange were obtained using 
total energy calculations for supercells with different spin 
arrangements, where the difference in the values of total 
energy originates only from spin degrees of freedom. To 
obtain the leading exchange integrals, we used a doubled 
cell with six Cu atoms. The supercell calculations were 
performed using the LSDA+ U method treating the cor- 
relation on a mean field level. This approach is necessary 
due to the Cu — O — Cu bond angle « 105° which leads 
to sizable FM contribution according to Goodenough- 
Kanamori- Anderson rules. 

The results are given in Tablc[I] {J^^^ was evaluated as 
the difference between the total J and J^^ = 4f^/{/e//)- 
The FM contributions significantly modify the size of the 
relevant exchange integrals, but preserve their AF nature. 
Here, we introduce the ratio a =Jd/Ji, which is zero in 
the simple kagome model and runs to infinity in case of 
decoupled chains. Certainly, a may depend on external 
parameters like the H position and the U values. The 
change in O — H distance drastically affects the Ji ex- 
change (Fig. [21 inset) in both compounds, and especially 
in haydeeite, where it becomes FM when the O — H bond 
is shorter than 0.8 A. Thus, further quantitative analysis 



is based on the empirical fact that total energy calcula- 
tions provide in general rather precise atomic positions. 
An accurate experimental determination of the H posi- 
tion is highly desirable for an improvement. The influ- 
ence of the Coulomb repulsion [J on a is much weaker and 
there are no drastic changes in verified region (Fig. d]): 
a is very close to 0.36 for kapellasite and stays in the 
vicinity of unity for haydeeite for the whole range of U 
studied. While GS's for a = and a = oo are relatively 
clear, the region in between is not studied. Therefore, we 
have performed exact diagonalization studies in order to 
clarify the influence a on the GS. 

It is well known that the classical GS of the pur e 
kagome HAFM (a = 0) is highly degenerate [1, HlISl- 
The additional diagonal bond Jd reduces this degener- 
acy drastically and selects non-coplanar GS's with twelve 
magnetic sublattices [s^j among the huge number of clas- 
sical kagome GS's. These classical GS's of the Ji-Jd 
model are characterized by a perfect antiparallel (Neel) 
spin alignment along the chains formed by diagonal 
bonds Jd and by a 120° spin arrangement on each tri- 
angle formed by NN bonds Ji. As a result, every two 
spin-sublattices are Neel-like antiparallel to each other 
and these two sublattices are perpendicular to one other 
group of two Neel-like sublattices. 

For the quantum model the GS and low-lying excita- 
tions have been calculated by Lanczos diagonalization for 
the finite lattice of iV = 36 considered previously in the 
literature for the pure kagome HAFM. Note that this fi- 
nite lattice fits to the magnetic structure of the classical 
GS. The calculated spin correlations (SqSr) for the clas- 
sical GS as well as for the quantum GS for a = 0.36 and 
a — 1.0 are shown in Fig. [5l For comparison we also 
show (SqSr) for the pure kagome system, i.e. Jd = 0. 
Obviously, the quantum GS spin correlation is drasti- 
cally changed by Jd- While for = the decay of 
the spin correlation function is extremely rapid, we find 
a well pronounced short-range order for a = 0.36 and 
a = 1.0 that corresponds to the classical magnetic struc- 
ture. This leads to the conclusion that even in the quan- 
tum model the GS has a non-coplanar magnetic struc- 
ture giving rise to enhanced chiral correlations. More- 
over, it is obvious from Fig. [5] that the magnetic corre- 
lations along the chains built by Jd bonds {R/Rnn = 2 
and 4) are strongest, indicating that the low-energy exci- 
tations might be S" = 1/2 spinons causing an effectively 
one-dimensional low-temperature physics similar to other 
2D models, e.g. the crossed-chain model [3^ and the 
anisotropic triangular lattice [s^. However, this issue 
needs further investigation. 

Finally we mention another important difference from 
the pure kagome system which is relevant for the low- 
temperature thermodynamics. For a = the singlet- 
triplet gap (spingap) is filled by 210 non-magnetic exci- 
leading to different low-temperature be- 
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havior of the specific heat C (power-law in T) and the 
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FIG. 5: (Color online) GS spin-spin correlation (SoSr) versus 
separation R= |Rj for the Ji-Jd HAFM on the kagome lattice. 
The results for the quantum 5 = 1/2 model are calculated for 
a finite lattice of A*' = 36 sites. The lines are guides for the 
eyes connecting data points. Note that for R = 2 and \/T2 
two non-equivalent spin-spin separations exist and the lines 
connect the points representing the stronger correlation. 



susceptibility x (exponential decay). By contrast we find 
that for a = 1 {a = 0.36) there are no (only a few) sin- 
glets within the spingap. Therefore we do not expect any 
basic difference in the low-T behavior of C and x- 

To summarize, we have performed electronic structure 
calculations for two new spin-1/2 kagome lattice com- 
pounds — kapellasite and haydeeite. Both compounds 
are 2D magnets, with two relevant AF exchanges: NN 
exchange Ji and the exchange along "diagonals" of a 
kagome lattice Jd- We find a =Jd/Ji ~ 0.36 for kapel- 
lasite and awl for haydeeite. The exchanges and thus 
a values are strongly dependent on the H position for 
which the experimental value is unlikely with respect to 
the total energy and should be reinvestigated. The pres- 
ence of significant Jd interaction leads to (i) non-coplanar 
magnetic order with 12 sublattices on the classical level 
which at least on a short-range scale, shows similarities 
to in the quantum model, and (ii) to the shift of the 
low-lying singlets out of the spingap. We especially em- 
phasize the crucial importance of Jd for all real materials 
with kagome geometry, which needs careful consideration 
in order to obtain the physically relevant model for the 
GS and the low lying excitations. This work is a starting 
point for study of these promising model compounds. 

Our predictions could be challenged and extended by 
low temperature experiments: neutron scattering and 
/iSR to probe the spin-spin correlation function, thermo- 
dynamic measurements (C and x, see above) to check for 
a possible spin gap, including pressure studies to modify 
a via a change of the 0-H distance. 

The work was supported by Emmy Noether program 
of DFG. We acknowledge fruitful discussions with S.- 
L. Drechsler and J. Schulcnburg. 
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